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Abstract 

The canonical quantisation of a charged particle in an intense laser background is developed 
and the associated two-point function calculated. It is shown how the causal propagator can 
be extracted. We provide an interpretation of this propagator in terms of degenerate quantum 
processes. The spectral properties of the propagator are derived and we show how multi-pole 
structures emerge as a perturbative feature. We also calculate the associated multiple wave- function 
renormalisations to all orders in perturbation theory. Our calculations are for both circular and 
linearly polarised plane wave laser backgrounds. 
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I. INTRODUCTION AND MOTIVATION 



Interest in the propagation of an electron in a laser background has a long history that 
significantly predates the recent advances in laser technology that is opening up the field 
to precise experimental verification. The exact solvability for the matter fields in a plane 
wave background was discovered by Volkov in 1935 [1], quite soon after the birth of quan- 
tum electrodynamics (QED). The advent of the first functional laser in 1960 [2] spurred 
extensive theoretical investigations into the effect of such a background on the propagation 



and interactions of matter 



Most notably, non-linear corrections to Compton scattering 
in such a background were calculated by various authors j^-^L, I n the mid 1990s experi- 



started to emerge which 



mental demonstration of such non-linear Compton scattering 
has revitalised theoretical interest in the field, see for example l9Hl3l. The commissioning 

n 

of several petawatt facilities worldwide, such as the Vulcan 10 petawatt project [14[ or the 
200 petawatt ambitions of ELI 15J, adds further impetus to the theoretical community to 
better understand quantum field theory in such intense backgrounds [if]]. 

Given the maturity of this field it is useful here to clarify our motivation for studying 
these problems afresh and to identify what new perspective we hope to bring to this area. 
Again, we need to look back to the early days of QED. Infra-red effects in QED have a similar 
antiquity to the above work on laser backgrounds, starting with the classic work of Bloch 



1* 



and Lee 



and Nordsieck [itI on soft infra-red effects followed by the work of Kinoshita 
and Nauenberg [19j on the collinear structures that arise both for massless charges and in the 
high intensity regime. The conclusion of these works can be summarised as follows: to arrive 
at a physical cross-section one must sum over all degenerate processes. This prescription 
for the infra-red underpins all attempts at calculating standard model cross-sections 20] 



but it has theoretical issues that are sti 
both initial and final state degeneracies 



1 hard to understand especially when there are 



2l|. In particular, the infra-red structures that 



arise from virtual corrections to the S-matrix need to cancel against a background of real 
degenerate (soft or collinear) processes that is indistinguishable from the original process. 
The existence of such degenerate backgrounds is, to say the least, problematic in most gauge 
theories, but a plane-wave laser provides a well defined background. So, coming from this 
infra-red perspective, a particle propagating through a laser is degenerate to one that emits 
and absorbs the same number of laser photons and thus these processes need to be summed 



over in order to describe the physical propagation of the charge. In this paper we want to 
develop this point of view and then back it up with a clear account of how such a system 
is quantised. Our hopes are that this will clarify some of the field theory underpinning this 
area of laser physics and, in turn, refine the techniques needed to tackle the infra-red sector 
of gauge theories, see also [22 1. 

After this introduction we will, in section UU develop a perturbative route to the Volkov 
propagator for a scalar particle propagating in a circularly polarised laser. This will be 
based on diagrammatic techniques that emerge from a systematic application of QED and 
the concept of degeneracy. Here we will calculate the leading and next to leading order 
contributions to the propagator. We will demonstrate the emergence of sideband modes 
within this perturbative framework already at lowest order in perturbation theory. Following 
that, in section II I II we shall take a completely different approach to the propagator. We 
will develop the equal-time quantisation of the Volkov field. This will allow for a precise 
identification of states and commutators and hence an all orders description of the Volkov 
propagator. We will see that these all orders results agree with the perturbative calculations. 
This will be followed in section HVl by a brief discussion of polarisation effects by now working 
in a linearly polarised laser background using both of the approaches taken earlier in the 
paper. Following the conclusion there are two appendices where some technical results are 
collected. 



II. SUMMING DEGENERATE PROCESSES 

In what follows we will consider perturbative corrections to the free propagator when a 
charged particle is in an intense laser background. We will not consider loop corrections but 
rather study the effects due to the absorption and emission of laser photons. In the spirit of 
the aforementioned Lee-Nauenberg theorem we will sum over all tree-level interactions with 
the laser which are indistinguishable from the propagator. We will not include in this paper 
any soft or collinear degeneracies. Such contributions would cancel against loop diagrams 
which we do not consider. This means that we will restrict ourselves to processes where 
the matter field has the same momentum, p, in both the initial and final states. Thus all 
emitted photons are necessarily degenerate with the laser and hence having four-momentum 
k. Furthermore, the number of emitted photons must be equal to the number of absorbed 
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laser photons as otherwise the matter field's initial and final momenta will not be the same. 
We will focus here initially on circularly polarised laser photons described by the potential 

A^(x) = a± cos(k-x) + sin(/c-x) , (1) 

which satisfies the gauge condition k-A = 0, where the null vector k will usually be taken 
to be pointing in the third direction: k^ = ko(l, 0, 0, 1). Furthermore, the polarisations 
vectors satisfy a\-a,2 = and a\ = a\ = a 2 < 0. For simplicity we will consider scalar matter 
interacting with the laser. The interaction Lagrangian has the form 

C int = -ieA» (V§^) + eMV . (2) 

To build up a perturbative description of the interaction with the laser background we use 
the Feynman diagram technique whereby we identify the basic interaction terms and their 
corresponding Feynman rules. For scalar matter two basic interactions arise corresponding 
to the two terms in the Lagrangian The cubic term corresponds to the expected three 
point interaction while the quartic term generates a four point vertex which is commonly 
called a seagull term. 

The four point vertex Feynman rule is obtained from the Fourier transform of the time- 
ordered product of free fields with the insertion of the four point interaction Lagrangian 
vertex 

J d 4 xe~ ip - x (0| T<P(x)^(0)ie 2 J d 4 za 2 ^(z)^(z) |0) (3) 
which after standard contractions generates 

fevf 2 * . Y- (4) 

The ie term in this expression arises through the use of time-ordering and encapsulates the 
causal nature of the propagator. At a more basic level, the significance of time ordering in 
(J3]) reflects basic commutators of the fields, which are well understood for this scalar field. 

The Feynman rule is then obtained from this by amputating the matter lines resulting in 
Fig. 1. Note that it is a special property of circularly polarised lasers, where A 2 = a 2 , that 
the only seagull term which survives is when one of the laser lines corresponds to absorbtion 
and one to emission and so both matter lines in (j3J) have momentum p. This is not the case 
for, e.g., linearly polarised lasers as we will see later in section HVl 
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FIG. 1. The seagull interaction 




= ie(ai + ia 2 )-p 




= ie{a\ — ia<i)-p 



p p + k 



p p — k 



(a) 



(b) 



FIG. 2. The three point interactions 



In a similar way, through insertion of the three point interaction Lagrangian vertex, one 
obtains the Feynman rules corresponding to Fig.'s 2a and b. 

We should here clarify the diagrammatic convention which we follow: laser lines coming 
from the left correspond to absorbtion of laser photons while lines going out to the right 
denote emission processes where the emitted photons are degenerate with the laser. 

The leading order corrections to the propagator are given by the seagull term of Fig. 1 and 
the diagrams in Fig. 3. We note that should the initial and final matter field have a different 
momentum, see, e.g., the diagrams of Fig. 4, then these contributions can appear in odd 
powers of the coupling e and it is clear that adding these processes to those of Fig. 3 would 
yield a result non-analytic in a. These are, anyhow, non-degenerate with the processes of 
Fig. 3 and cannot contribute to the propagator. 

In this spirit we now consider the leading and next to leading perturbative corrections to 
the free propagator with the matter field having the same initial and final momentum. 




p p + k p 



p p - k p 



(a) 



(b) 



FIG. 3. Leading order contributions to the propagator 
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p p+k 



p p p+k 



p p + k p + 2k p+k 



(a) (b) (c) 

FIG. 4. Some leading and next to leading non-degenerate laser-particle interactions 

A. The propagator at order e 2 

The leading order corrections to the free propagator are given by Fig. 1 and the two 
diagrams of Fig. 3. Fig. 1, with the four-point vertex, contributes 

^-r 1 —, ■ (5) 



Note that in this, and for the rest of this section, we have suppressed the it term for 
compactness of presentation. Expression ([5]) is a double pole in p 2 — m 2 and this corresponds 
to a mass shift effect. Writing the propagator as the free propagator plus tree diagram 
corrections 



p2 _ m 2 _ fi m 2 p2 _ m 2 
we can immediately identify the mass shift as 



1 + -5 = (-i6mf) + ... 

pA _ m A 



(6) 



Sm 2 = - e 2 a 2 , (7) 
which is a positive correction as a 2 is negative. This is the mass shift that was first identified 



by Sengupta 



23| and later rediscovered by various authors , that arises when a charge 



is in such a laser background resulting in an effective mass m* where 

ml = m 2 + e 2 (-a 2 ). (8) 

The matter lines between the absorbtion and emission vertices in the two diagrams in 
Fig. 3 carry momentum p + k. Putting the diagrams together we have 

_ g 2 f i \ (! av p\ 2 -f (a 2 -p) 2 ) ( - I - 

\p 2 — m? J \{p + k) 2 — m? (p — k) 2 — m? 

For convenience we define 



2 (avP) 2 + (arpf , 1fV . 
W ? = (fW ' (10) 



and note that when the terms in the large brackets in ([9]) are placed on a common denomi- 
nator the resulting numerator cancels one of the poles in p 2 — m 2 yielding the contribution 

2ie 2 w 2 ( p ' fc ) 2 . (11) 

p {p 2 — m 2 ) ((p + k) 2 — m 2 ) ({p — k) 2 — m 2 ) 

At this stage we use partial fractions to disentangle the fraction above: 

{p-kf 1 

(p 2 — m 2 ) ((p + k) 2 — m 2 ) ((p — k) 2 — m 2 ) 4(p 2 — m 2 ) 

1 1 

+ 7777 TTT; 77 + 



]((p+k) 2 -m 2 ) 8((p - k) 2 - m 2 ) ' 
We identify three structures here. The first is a single pole in p 2 — m 2 and thus a wave- 
function renormalisation effect on the initial pole. The other two terms are new poles at 
{p ± fc) 2 — m 2 . These poles first appear at order e 2 and are of exactly the form obtained 



m 
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26] which were obtained using non-perturbative arguments. We emphasise that they 
are completely perturbative and that they appear here from the lowest order correction to 
the free propagator. We will see in a moment that they also acquire the mass shift of (JTj) 
from diagrams entering at order e 4 . We will write the terms involving poles in p 2 — m 2 in 
the form 

,■7(0) „- 7 (0) r • i 



1 + 9 - „ (-i5m 2 ) + ... 



p 2 — m 2 — 5m 2 p 2 — m 2 
where the bracketed superscript in indicates that this corresponds to a denominator 
p 2 — m 2 = (p + Ok) 2 — m 2 . From this equation we see that at this order the residue of the 
pole is 

zfK) = l-e 2 ^, (13) 

which is a finite renormalisation. 

The new pole structures at {p ± nk) 2 — m 2 first appear, for integer n, at order e 2n and 
the appropriate wave function renormalisation type constant, Z% n \ will be a power series 
in the coupling starting at order e 2n . From the above we have at leading order 
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There is an important point to note here that we will return to later in the paper. In 
QED without a background the wave-function renormalisation constant Zi is momentum 
independent. This is because it is identified with the residue of the appropriate pole, i.e, it is 
evaluated on-shell. Relativistic invariance then dictates that it could only be a function of p 2 



Zt\w p ) = e 2 ^ . (14) 



which becomes m? in Z2 irrespective of the choice of positive or negative energy pole. With 
a background things are different as we have an additional vector k that can, and indeed 
does, enter into the calculations. This means that the wave-function renormalisations are 
now not constants but can be functions of p-k, p-a\ and p-a-i, with p on-shell. In addition, 
the residue will be different depending on whether the positive or negative energy pole is 
picked up. 



B. The propagator at order e 



We now proceed to the next order correction. These are given by the thirteen diagrams 
of Fig. EJ 

We note from these diagrams that the internal matter lines may carry momentum p, p±k 
and, now, p ± 2k. Thus we can expect the poles from the lowest order diagrams to receive 
corrections and additionally new poles at (p ± 2k) 2 — m 2 to arise. It is helpful to divide 
up the study of these diagrams according to the number of four point and of three point 
vertices in the diagrams. Diagram Fig. [5^ which only contains the four point vertex only 
features poles in p 2 — m 2 . It simply generates the next order term in the geometric series 
of the mass shift found above in ([6]). In other words ([7]) is not altered by a higher order 
correction: 
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p2 _ m 2 _ fi m 2 p2 _ m 2 



p2 _ m 2 



-i5m 2 



p2 _ m 2 



:-iSm 2 ) 2 



■ (15) 



The next diagrams, Fig. |5b-g, with a mix of four point and three point vertices only 
generate double poles (some single poles cancel in the sum of these diagrams). They yield 



^a 2 w 2 v 



p-k 



.p-k 



p2 _ m 2 



+ i 



p-k 



+ i 



p-k 



77T 



• (16) 



(p + k) 2 — m 2 J 4 \(p — k) 2 

Finally at this order, the remaining diagrams of Fig. [51 without four point vertices, only gen- 
erate single pole structures (with double pole terms cancelling in the sum of these diagrams). 
They yield 

3 i 1 i 1 i 



4 4 
e w p 



32 p 2 — m 2 16 (p + k) 2 — m 2 16 (p — k) 2 — m 2 

1 % 

+ 



+ 



1 



64 (p + 2k) 2 - m 2 64 (p - 2k) 2 - m 
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p p p 



p p p — k p 



p p — k p p 



(a) 



(b) 



(c) 




p p p + k p 
(d) 





p p — k p — k p 



(g) 




-> — * >- 



-> * — >- 



p p + k p + 2k p + k p 



(k) 





p p + k p p + k p 
(h) 




p p+k p p—k p p p—k p p + k p 

(i) 0) 





-> — * — ► 



p p—k p p—k p 



(1) 




p p — kp — lkp — k p 



FIG. 5. The 0{e i ) degenerate contributions to the propagator 



These last two equations clearly display the behaviour stated in Reiss 25] that the four 
point vertices solely generate the mass shift while the three point vertices yield wave function 
renormalisation structures. We will also see how this behaviour arises in the next section 
from a canonical quantisation approach. However, later in section IIVt we will see that this 



divide does not hold for linearly polarised lasers. 

The cross terms in (Tl2|) involving the order e 2 mass shift and Z^ 1 ^ can be recognised in 
f fT6|) showing again that the mass shift is determined by the lowest order terms and receives 
no additional corrections. In the next section we will understand at all orders why these 
new poles acquire the same mass shift ([7]). The remaining terms above tell us that 

that the next correction to the (p ± k) 2 — m 2 wave function renormalisation is 

2 4 

^(±1)/ \ 2 W P A W P / 1Q \ 

Z\ \w p ) = e — -e — , (18) 
and that the first term of the wave function renormalisation (p ± 2k) 2 — m 2 is 

4 ±2> K)=e 4 g. (19) 

We will proceed to obtain a better understanding of all of these structures in the next section 
when we investigate the two-point function for the Volkov field. 



III. TIME ORDERED PRODUCTS AND THE VOLKOV PROPAGATOR 

In this section we want to analyse in a more fundamental way the Volkov propagator 
and, in so doing, we will derive the all-orders generalisation of the perturbative results 
obtained in the previous section using the Feynman diagram method. Again, we consider a 
circularly polarised (plane-wave) laser described by the classical potential (OQ) and identify 
the (complex) Volkov quantum field <f) v ( X^j clS db solution to the background laser coupled 
Klein-Gordon equation: 

(D 2 + m 2 )4 >v (x)=0, (20) 

with D 2 = D^D^ and our conventions are such that = + ieA^. 

In order to understand how (120]) can be solved it is useful to recall the method used in 
the free case where we wish to solve the Klein-Gordon equation (d 2 + m 2 )0 = with now 
d 2 = d^d^. The solution to this is built up by noting that the plane waves e ±ip ' x satisfy 
the equation Q 2 q ±1 p x = _p2 e ±jpx g Q ^ ^ four momentum p is taken to be on-shell: 
p — >■ p = (E p ,p) where E p = ^j\p\ 2 + m 2 so that p 2 = m 2 , then a superposition of these 
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plane-waves gives us the general solution to the Klein-Gordon equation 

f3 



/^(e-^a(p)+e^(p)) , (21) 

where we have adopted the convenient notation that d 3 p = d 3 p/(2n) 3 . This is a quantum 
field because the modes a(p) and tf(p) are quantum operators whose commutator properties 
are central to the construction of the corresponding states and their propagation. 

In the Volkov case a similar procedure is possible. That is, we can replace, for example, 
q-w-x by a function E(x,p) where D 2 E(x,p) = —p 2 E(x,p) (see, for example, 27]). So again, 
if the momentum is taken on-shell as above, then we can build up a solution to (1201) using 
E(x,p). But, in the light of our discussions in the previous section, we expect that the 
on-shell condition in such a laser background should involve the effective mass (IE]). This 
suggests that we should really be looking for a replacement of the plane wave by a distorted 
field D(x,p) such that now 

D 2 D(x, p) = (-J9 2 + 5m 2 )D{x, p) . (22) 

Then if we go on-shell at the starred- mass: p — >■ p with now p = (E*,p) where E* = 
\p\ 2 + ml so that p 2 = ml = m 2 + 5m 2 , we will still satisfy the coupled Klein-Gordon 



equation with mass m: 

D 2 D(x,p) = (-p 2 + 5m 2 )D(x,p) = -m 2 D(x,p) , (23) 

and in this way build up the Volkov field. Obviously these two approaches are not unrelated 
and indeed we have 

/ Sm 2 \ 

E(x,p) = exp —i -k-x D(x,p) . (24) 



2p-k 

But, as we shall see, the transition from on-shell to off-shell momentum is a subtle point in 
our analysis and focussing on the distortion field D(x,p) will make this step much clearer. 

Specialising now to the circular polarisation case, where we write the distortion field as 
D c (x,p), we find that 

D {x p) — e~ ip ' x e ie( - wi p sin ( k - x )+ w 2p cos(k-x)) ^25) 

where 

p-ai p-a 2 . . 

Wlp - = -— and w 2p =-^. (26) 
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Using these, the equation of motion (1201) can then be exactly solved and we find that the 
Volkov field can be written as 

M*) = (D c (x,p)a v (p) + D c (x, • (27) 

The momentum assignment in (1271) might seem a little unnatural as we could equally well 
satisfy the equation of motion f )20|) by using D c (x,pytf y (p) in the second term. But, as we 
shall see, this would then obstruct the gauge transformation properties of the Volkov field. 

To understand this it is useful to note that there is a residual gauge freedom in our 
description of the background vector potential and we could equally well use A^x) + 
8^9 (x) where 

9{x) = 0i sin(fc-x) + 6 2 cos(k-x) , (28) 

with 6*i and 6 2 arbitrary constants. This gauge transformation can equivalently be written 
in terms of the polarisation vectors as 

+ and a£->a£-0 2 F- (29) 

Then, the terms (1261) have particularly simple gauge transformation properties and under 
(129|) we have 

Wip — > Wip — Q\ while w 2p — > w 2p — 2 . (30) 

Hence, from (1251) and fl30l) . we see that the Volkov field <p v (x) as defined in (1271) will transform 
under a background field gauge transformation as expected for such a matter term: 

^ v (x)^e-^V v (a:). (31) 

As in the free theory, we would like to interpret the Volkov field's modes that enter into 
(1271) as creation and annihilation operators and hence develop a description of the Volkov 
states and thence the Volkov propagator. In order to do this we need to be able to invert (|27|) 
so that we can express the operators a v (p), etc, in terms of the Volkov fields. In the free case, 
this inversion is achieved using basic Fourier transform ideas, but now the distorted plane 
waves complicate things. However, simple physical guidelines will allow us to generalise the 
Fourier inversion to this Volkov description. 

The guiding principles in finding an inversion for (|27|) are that the Volkov modes are 
gauge invariant and time- independent. Given that the inversion result should reduce to the 
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free field expression when the coupling is switched off, after a little experimentation we are 
lead to the expression: 

a v (p) = -i J d 3 x ( (D D c (x,p)) f <p v (x) - D c (x,p) j D <f) v (x)) , (32) 

with similar expressions for the other modes. In this equation gauge invariance is ensured 
through the use of the covariant derivative D = d + ieA and time independence follows 
from the equations of motion ( 1201) and ( 1231) along with the simple identity that OqDq = 
D 2 — DiD 1 — ieA D . The actual proof of (1321) . and its generalisation to the other modes 
and the associated creation operators follows from the identity that 

-i J d'x ((D D c (x,p))^ v (x) - D c (x,pyD <P v (x)) (33) 

= j% {mq)a v (q) + I(p,-q)bi(q)) , 

where we have defined the integral 

I(p,q) := -i J d 3 x ((D D c (x,p)) j D c (x,q) - D c (x,pyD D c (x,q)~) . (34) 

In Appendix [B] we show that this can be evaluated explicitly with the result that 

I(p, q) = (2tt) 3 (E; + E* q )5\p - q) . (35) 

The proof of the inversion result (|32|) then follows immediately. 
Given the expression ( |32l) for the Volkov modes and its adjoint 

4(g) = i J d 3 y ((AAMKfo) - D c (y,q)(D cP v (y)y) , (36) 

we can exploit the time-independence of these results to calculate their commutator by using 
the equal-time canonical commutators of the Volkov fields. From this we see that 

[a v (p),<(<l)] = (^) 3 2E;5 3 (p-q). (37) 

Given these commutators we can now construct the unequal-time commutators: 

fev(*),4(y)] = j d ^{D c (x,p)D c (y,py - D a (x,-p)D e (y,-tf) , (38) 

the study of which lies at the heart of the causal properties of these fields and will be 
the subject of another publication. For now we shall consider the Volkov states and their 
propagation. 
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From (|37|) we see that the Volkov modes satisfy the usual commutation relations which 
allows us to identify them as creation/annihilation operators, but with a relativistic nor- 
malisation given by the effective energy E* So acting on the Volkov vacuum, where 
a v(p)|0)v = ^v(p)|0)v = 0, we can construct multi-particle Volkov states. However, in 
contrast to the free case where 0'(x)|O) can be understood as a superposition of single par- 
ticle states since (d 2 + m 2 )(j)\x)\0) = 0, for the Volkov field 0^ (x)|0) v is a much richer 
ensemble of states. As discussed in Appendix A, the distortion term can be written as a 
superposition of Bessel functions and thus we can write 



oo 



4>U x )\°)v = E *l(*)l°>v ( 39 ) 

where 



n=— oo 



4(*)|0)v = /||^e t(p "-" fc) -e-^J n (e Wp -)a v t( P )|0) v . (40) 
These now satisfy a shifted free equation of motion: 

(( 5 + mA; ) 2 + m 2)0l i ( x )|o) v = O. (41) 



We shall refer to them as sideband states in analogy with Reiss' 28] discussion of the multi- 
pole structure of the Volkov propagator. 

As is the case for the free theory, we expect the causal propagator for the Volkov field to 
be related to the two-point function 

v (O|T0 v (x)4(y)|O) v , (42) 

where the time ordering is given by 

T<f> v (x)<f^(y) = 0(x° - y°)0v(aO4(y) + %° ~ *°)4(v)0v(*) (43) 
and 6{z) is the step function which we can write as 

poo ■ 

9{z) = - lim / dX- e iXz . (44) 

J_ 00 A — ze 

The mode decomposition fl39l) of the Volkov field will induce a double decomposition for the 
two-point function: 

v <O|T0 v (a:)0i(y)|O> v = Ev(O|T0 v (^)0i(2/)|O)7 := £ v (O|T0 m (x)0t(y)|O> v (45) 



m,n 
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and this poses a challenge to the direct identification of the propagator. 

We recall that in the free theory the field acting on the vacuum creates particle states 
which can then be annihilated at a later time by the conjugate field. This insight allows us 
to identify the propagator as (the vacuum expectation value of) the time-ordered product 
of these fields. For the Volkov field we have seen that this particle interpretation is very 
different. The modes of the field, </4(x), can be interpreted as creating particle states (1401) . 
However, each of these modes has its own on-shell type of condition fHTj) . Thus the annihila- 
tion of such particles at a later time requires the corresponding conjugate mode. Therefore, 
in the sum (145 p . only the diagonal terms have an interpretation as a propagator. This iden- 
tification is the operator version of the degeneracy argument given in the previous section. 
The interpretation of the off-diagonal terms in the two-point function is less clear but is 
related to the diagrams in Fig. 4 and will be addressed elsewhere. 

From the commutation relations ( 137|) and the expression (T40|) we thus identify the prop- 
agator as a sum over terms of the form 

v (O|T0 n (x)0t(i/)|O) v = - j d ^J\^^{e^- E i +nk ^°-y°) e^--^-y) (46) 



+ e - i{ X-E* p -nk )(x°-y°) e -i(p+nk).(x-y)\ j2 ^ 



In order to make this look more like a propagator we need to go from on-shell expressions 
for the momenta to off-shell ones. The way this is done is intuitively clear but the details 
are a little subtle. Basically we want to identify the integration variable A as the missing 
four-momentum component p°. To see how this works in practice, we note that this integral 
is a sum of two integrals and we need to make different changes of variables in each. In the 
first integral we set p° = A — E* + nk° while in the second — p° = A — E * — nk° to give: 



MTM*WM\0)v = -} 3 4 pe ip ° {x °- y0) ( p0 - E; z nk o _ l£ e*fe-*-M (47) 



+ 



E;-p° + nk° -ie " J 2E*" n 



We now shift the 3-momenta, and again we have a different shift in each term: in the first 
integral we change p — > —p + nk while in the second p — » p — nk. 



In both cases (since E* = E*_ p ) we get the shift E* ->■ E*_ nk = J\p - nk\ 2 + mj, the 
change to Wp are, however, a little more subtle as we now explain. Under our changes of 
variable the on-shell momenta change in different ways. In particular, under the first shift 
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we get p — > (E*_ p+nk , —p + nk) while under the second we have p —> (E*_ nk ,p — nk). Both 
of these are on-shell but at different momenta. However, exploiting the fact that Wp = W-p, 



we can view these shifts within Wp as 



p^p-nk = (-E*_ nk , p-nk) and p p - nk = (E*_ nk , p-nk). (48) 

These are the negative and positive on-shell momenta corresponding to the on-shell condition 
(p — nk) 2 = m 2 . 

Under these changes of variables, ( 1471) becomes (after putting terms over a common 
denominator) 

V <0 W „ W «M|0> V =/ ^ (p _ nfc)2 '_ m; + i6 e^» (49) 



OP* ^ p-nk J ^ OP* °n\^ UJ p-nk, 

LCj p-nk LCj p-nk 

We now note that the terms in brackets greatly simplify when evaluated at the two poles of 
the overall denominator where they become either ^(ew-^.) or J^( eu p^fc)- The off-shell 
version of these are both then J^(ew p _ n fc) = J^(ew p ) since a-k = k 2 = 0. This allows us to 
write the full Volkov propagator in a fully off-shell form as (noting that J 2 (a) = J 2 n (a)) 

V [ d 4 p l ^ ew p) c iv(x-y) _ s 5Q \ 
J (p + nk) 2 — m 2 + ie 

n=-oo u xr ' * 

From this expression we can deduce that the all-orders expression for the laser induced wave 
function renormalisation is 

Zi n \w p )=3 2 n (ew p ). (51) 

Using this identification and the standard expansions for the Bessel function, we see that to 
order e 6 







= 1 e 2W he^< 
2 32 


— e 




[w p ) 


2 4 
W HI 

= e — — e — + e 
4 16 


768 




[w P ) 


4 6 
W 111 

- e 4 p e 6 p 
64 384 ' 






[w p ) 


w 6 

- e 6 p 
2304 ' 





.5w p 



6 

576 ' 

(52) 



Here we see confirmation and an extension of the perturbative calculations presented in 
section 2. Indeed, the identification (I5T|) allows us to conclude that to all orders, 

Zf\w v ) = Zt n \w p ) (53) 
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•2 



p V p p + 2k p p — 2k 

FIG. 6. The four point Feynman rules for linear polarisation 

= tea ■ p 



p p + k p p — k 

FIG. 7. The three point Feynman rules for linear polarisation 

IV. POLARISATION DEPENDENCE 

In this section we will consider the scalar matter in a linearly polarised laser background 
to explore the polarisation dependence of our results developed in the previous two sections. 
For a linearly polarised laser the background potential becomes 

A^(x) = a^cos(k-x) , (54) 

where the amplitude vector satisfies k-a = 0. Following the techniques of Sect. II, we find 
the Feynman rules shown in Fig.'s[6]and[71 

Comparing with Fig. 's [Hand [2] we see ^ na * the Feynman rules are polarisation dependent. 
The most important difference is the existence of additional four point vertices where both 
photons are incoming or both outgoing. At order e 2 these rules do not generate diagrams 
in the propagator as the initial and final matter momenta would not be the same. Thus the 
propagator at this order is still given by the degenerate diagrams of Fig. 's Q] and El Thus 
the mass shift is now given by 



5m 2 = -^~, (55) 



where the factor of a half follows from the difference between the Feynman rules in this 
case. Here we see the polarisation dependence of the mass shift which is an example of the 



dependence of the mass on the form of the laser field 29] . 
The wave function renormalisations similarly are 

Z (0) =1 _ e 2^ (5g) 

where 

P ' a (57) 



p ■ k 
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V p + 2k P P p + 2k p + k P 

FIG. 8. Examples of additional, higher order, degenerate diagrams in a linearly polarised laser 



and 



2^ = *\. (58) 



At order e 4 there are twenty one degenerate diagrams that contribute to the propagator. 
These are those of Fig. plus an additional eight diagrams which involve the new seagull 
Feynman rules. These include, for example, the diagrams of Fig. [SJ 

We will simply summarise the results obtained for the propagator up to order e 4 . As 
with the circular polarisation case, we find that the mass shift is not altered by higher order 
diagrams. The result is still that of ( 15 5 p . However, for the wave function renormalisations, 
with notation being a simple extension of the circular polarisation case, the results are rather 
different. Defining 



a 2 



v P = 5— r , (59) 
8p ■ k 



we obtain 



and 



Z f( v „) = l-e*f-e'(f-g), (60) 



while 



4 ±1 W^)=^-e'(^±^), (61) 



^ ±2 W%) = e 4 (! + f±^l- (62) 



Again we note that is of order e 2n and higher. Furthermore we see that although once 

again only the seagull vertex produces contributions to the mass shift, the wave function 
renormalisation for linear polarisations includes contributions from both seagull and three 
point vertices. Thus the fact that the wave function renormalisation for circular polarised 
lasers was generated by three point vertices alone is demonstrated to be a polarisation 
artefact. We also note that in contrast to the circular polarisation case, Z% does not equal 



Z 2 



(-n) 
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At the operator level, the distortion term for a linearly polarised background is 

D p) = Q- i P- x e K e ^ ai - n (.k-x)+e?vp8in(2k-x)) 

oo 

= e -i P -x e ink - x J n {eu p ,e\). (63) 

n=— oo 

An important point to note about this expression, and one that did not arise for circular 
polarisation, is that the arguments of the generalised Bessel function are not invariant under 
the replacement p — > —p. In particular, although U-p = Up, we now have from ( 159|) that 
V-p = —Vp. This means that in the derivation of the propagator for a linearly polarised 
background, we need to replace equation f l4"6"j) by 

v (o\T<t> n (x)<t>UyMv = - j d ^^j^( el(X ' E;+nk0){x °~ v0) eJ( -" n - >fe "- )j n(^ P -,e 2 ^) 

(64) 

+ e - l{ X-E;-nk )( X °-y°) e -^nk)i,-yjj2^ eu ^_ e 2 v _^ 

If we now apply the same arguments as before for circular polarisation, we obtain in place 
of ( H9|) the expression 

v ( Ej-nk ~ (P° ~ l2 , _ 2 , , E* p _ nk + {jP-nk°) i2f 2 x 

X l OP* J n[ eU f^j,, eV f^ k )+ 9F * J n{ eU p-nk' e V p-nk) ) • 

V Zrj p-nk LCj p-nk 7 

Repeating the argument under (H9|) and using the identity (lAllj) . we obtain the propagator 
in a linearly polarised background 

y / iJ "7 e y e^~^. (66) 

■ / -^' / (» + nkY — mi + ie 

n=-oo ^ v ^ ' * 

Hence the wave function renormalisation constants in this linearly polarised background for 
all the states at all orders are given by 

Zf\eu p , e 2 v p ) = Jl(eu p , e 2 v p ) . (67) 

It is easy to verify that this agrees with the perturbative results (160!) to (|62|) . 



V. CONCLUSIONS 



In this paper we have calculated the propagator of an electron in a laser background in two 
different ways as well as in different laser backgrounds. For simplicity this work was carried 
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out in scalar QED. The first approach we pursued was motivated by infra-red physics: the 
physical propagator is composed of the sum of all degenerate (indistinguishable) processes. 
For example, an electron which passes freely through a laser is indistinguishable from one 
which, say, absorbs a laser photon and emits a photon of the same energy and momentum 
back into the laser. Thus the propagator in such a background should involve a sum of tree 
diagrams with Fig. 's[T] and |3] showing the order e 2 degenerate processes. 

We have thus shown that this leads to the well known mass shift and also that the 
perturbative corrections generate additional poles (sideband states) in the propagator of the 
form 

1 (68) 



(p ± nk) 2 — ml ' I — | 
for integer n. We stress that sideband states arise perturbatively [28], some already from the 

lowest order corrections to the free propagator. The wave function renormalisation factors, 

i.e., the numerators of these sideband poles, start at order e 2n (with n as above). In this way 

we have calculated the degenerate corrections to the free propagator up to next to leading 

order. 

In the scalar theory we consider here, there are both three and four point interactions 
between the matter and laser fields. The form of these Feynman rules depends upon the laser 
polarisation. However, for both circularly and linearly polarised lasers the mass shift was 
solely generated by contributions from the four-point vertex at lowest order. All higher order 
diagrams which generate double poles simply yield the higher order terms in the geometric 
expansion (|T2|) . It was also shown that the mass shift in all poles was the same. 

The wave function renormalisations for circularly polarised lasers only involve the three 
point vertex structures. There is also the symmetry that Z^' — Z% n \ However, for linearly 
polarised lasers, the wave function renormalisation structures depended upon both sets of 
vertices. Another new feature of the linearly polarised case is that 4 n) ^ Z ( f n \ i.e., the 
symmetry between these sideband states visible for circularly polarised laser backgrounds is 
removed for linear polarisation. 

The second approach which we have pursued in this paper was based on the canonical 
quantisation of these fields. We are not aware of such an investigation in the literature 



although we note |30[ and some lightcone quantisation work [3jJ, l32[- Our approach yielded 
an all orders result which completely agreed with the Feynman diagram results. We have 
constructed the canonical quantisation for Volkov states including their unequal time com- 
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mutators. The two-point function for Volkov fields was written as a double sum over creation 
and annihilation type modes. These modes satisfy a free equation of motion, see (|4ip . but 
shifted in accord with the sideband states discussed above. However, as only the diagonal 
terms in the sum ( I45p have a propagator interpretation we restricted our attention to these 
terms. The Volkov propagator defined in this way could be expressed as a sum over sideband 
states and the normalisation given in terms of Bessel functions (generalised Bessel functions 
for linearly polarised laser backgrounds) see fl5Tj) and fl67|) . Expanding these results yielded 
all of our previous perturbative calculations. It is not immediately obvious to us that our 
infra-red motivated sum over degenerate diagrams and the diagonal terms in the canonical 
quantisation double sum had to agree like this. 

It is useful to recall here that, even without a background field, the pole in the free prop- 
agator is modified in perturbative calculations to yield the full spectral decomposition of 
the two-point function. This now includes renormalised poles, bound states and branch cuts 
corresponding to scattering states. What we have shown in this paper, without calculating 
any loop corrections, is that interactions with the background induce a modified spectral 
structure with infinitely many discrete poles corresponding to the sideband states. It is in- 
triguing to speculate how the loop corrections will further modify these sideband structures. 
In particular, one might ask whether infra-red divergences will prevent sideband states from 
displaying mass gaps and, following Lee-Nauenberg, what real processes should be summed 
over to construct finite inclusive cross-sections in a laser background. 

This work raises many further questions such as the extension to fermionic QED, the role 
of the off-diagonal two-point function structures and the general form of the vertices both 
perturbatively and non-perturbatively. They will be addressed in future publications. 

Appendix A: Generalisations of Bessel functions 

The distortion of the plane wave solutions due to the interaction with the laser are best 
described using Bessel functions. For a circularly polarised laser there are two ways to 
formulate this. One is in terms of standard Bessel functions while in the other a generalised 
formulation of Bessel functions is introduced. It is useful to use both of these and understand 
how to move between them. For a linearly polarised laser the only option is the generalised 
route and this will be easier to understand once we have discussed the circular case. 
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1. Circular Polarisation 



In the expression (1251) for the distortion term D c (x,p) we can write 

wipsm(k-x) + W2pCos(k-x) = Wp sxn(k-x + i/jp) (Al) 

where 

Q Q Q IjU^ / A \ 

w ip + w 2p = w p an d tanyjp = . (A2) 

Using this we can then write the distortion term in terms of standard Bessel functions: 

00 

D c (x,p)=e-^ x e m{k ' x+ ^ p) Jn(ewp). (A3) 

n=— oo 

The advantage of this formulation is that in the propagator we will always have even powers 
of and the phase ipp will cancel so that any ambiguity in the sign of can be avoided. 
But, in deriving the inversion formula (I3l)j) . odd powers of will arise and it is much 
simpler to use a generalisation of the usual Bessel functions as outlined below. 

Rather than rewriting the sine and cosine terms in D c (x,p) as a shifted sine as in OA II) 
and then using the Bessel representation, we can directly express both the sine and cosine 
terms as Bessel functions to arrive at the representation 



00 



D c (x,p) = e"^ e mk x Uew lp} ew 2p ) , (A4) 

n=— 00 

where we have defined the generalised Bessel function J n (a,(3) by 

00 

J„(a,/3)= * r Jn-r(a)J r (/3) • (A5) 

r=— 00 

These Bessel functions satisfy various identities that follow quite easily from their definition 
and we note, in particular, that J n (0,0) = 5 n o and, most importantly 

2nJ n (a,(3) = a(J n _i(a, 0) + J„+i(a, 0)) + i/3(J n _i(a, 0) -J n+1 (a,/3)) . (A6) 

As we shall see, both (1A3I) and ( 1A4I) are useful representations of the distortion factor and 
we can shift between them by using the identity 

e in ^J n (ew p ) = J n (ew lp , ew 2p ) . (A7) 
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2. Linear Polarisation 



In the linear case the only option is to expand both the sin(fc-x) and the sin(2fc-x) terms 
in D e (x,p) as series of Bessel functions to give the generalised Bessel function expansion 
(see Appendix J of [33j) 



oo 



D t (x,p) = e"** e mk ' X Ueu P , e 2 v p ) , (A8) 

n=~ oo 

where now the generalised Bessel functions J n (a,(3) are defined by 

oo 

J„(«,/3)= ^ J n _ 2r (o)J r (/3) . (A9) 

r=— oo 

Again, these Bessel functions also satisfy J n (0, 0) = 5 n Q and now 

2nJ n {a, /3) = a(J„_i(a, /3) + J„+i(a, /?)) + 2/3(J n _ 1 (a:, /3) + J„+i(a, /?)) . (A10) 
We also note that the generalised Bessel functions satisfy the identity 

J n («,-/3) = (-l)V_ n («,/3). (All) 

Appendix B: Volkov Inversion Formula 

Here we shall derive (|35|) for circular polarised Volkov fields, the linear polarised result 



follows in much the same way (see ( 34|) for the full details) and note that, for simplicity, 
in this section we are setting the coupling to unity 

The first observation to make is that the gauge invariance of can be exploited to 
work in the gauge where, in addition to k-A = 0, we have A = A 3 = 0. This means that 
the covariant derivatives in (1341) can be replace by ordinary derivatives. Using (1251) it is then 
straightforward to see that 

(Bl) 

Writing the final product in terms of the generalised Bessel functions (1A5|) 



D c (x,pyD c (x,q) = e^ x ^ e m ^J„K 9 - - w lp , w 2q - w 2p ) , (B2) 
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and using the standard exponential description of the sine and cosine terms which allows us 
to perform the x-integral, we arrive at the expression 

oo 

I(p, q) = (2tt) 3 S 3 (p-q + nk)e i{E ^ +nko)x ° 

n=—oo 

*((E* + E*)H n - ±k (w lp +W 1 q)(j n - 1 + J n +l) - %k (w 2 p + W 2 q)Qn-l - Jn+l)J , 

where the arguments of the Bessel functions have been suppressed but are all the same 
differences as given in (1B2|) . We now need to consider separately the cases n = and n ^ 0. 

When n = the delta function implies that E* = E* and hence Wi p = Wig as well as 
w 2 p = W2q, which means that the arguments of the Bessel functions vanish. Now using the 
result that J n (0, 0) = 6 n o we get the contribution to I(p, q) of the term (27r) 3 5 3 (p— q)(E*+E*) 
which is the right hand side of the identity ( |35|) that we wish to show. Hence we now need 
to verify that when n^O we get no further contributions to I(p, q). 

When n / the delta function and gauge condition imply that both w\p ^ W\q and 
W2p 7^ so we can write the sum as 

6 3 (p -q + nk)e iiE r- E « +nko)x ° 

x Ue; + E* q )\ - \k Q ( w W + w iA ( w _ w^X^ + J n+1 ) (B3) 
v \wiq-w lp J 

-\k ( W2P + W29 ) i{w 2 q ~ W 2 p)(j n -l "Jn+l)) • 
\W 2 q-W 2 pJ / 

The delta function in this also implies that both p-a\ = q-a\ and p-a 2 = q-a 2 . Using this we 
see that 

Wip + w li! _ w 2p + w 2 q _ q-k + p-k _ E* + E* - 2p 3 - nk 
wiq-wip w 2 q-w 2p p-k - q-k E*-E* + nk 

Using this and the identity (IA6j) we can write the sum (1B3|) as 

£ S\ - j + nk)e«°--^*° (( E ; + El) - \K ( 3 ^l'^ ~ k ) *•) *. ■ (») 

Putting the bracketed terms over a common denominator quickly shows that this sum van- 
ishes as required. 
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